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Abstract 

The exact global solution of the Einstein equations [Neugebauer & 
Meinel, Phys. Rev. Lett. 75 (1995) 3046] describing a rigidly rotating, 
self-gravitating disk is discussed. The underlying matter model is a per- 
fect fluid in the limit of vanishing pressure. The solution represents the 
general-relativistic analogue of the classical Maclaurin disk. It was derived 
by applying solution techniques from soliton theory to the axisymmetric, 
stationary vacuum Einstein equations. In contrast to the Newtonian so- 
lution, there exists an upper limit for the total mass of the disk - if the 
angular momentum is fixed. At this limit, a transition to a rotating black 
hole, i.e., to the Kerr solution occurs. Another limiting procedure leads 
to an interesting cosmological solution. These results prove conjectures 
formulated by Bardeen and Wagoner more than twenty-five years ago. 
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1 Introduction: 

Rotating bodies in general relativity 

The problem of self-gravitating rotating bodies is a global one: One has to solve 
simultaneously the interior equations (with matter) as well as the exterior (vac- 
uum) equations. The shape of the surface cannot be prescibed arbitrarily. A 
famous solution within Newton's theory of gravitation is the Maclaurin spheroid 
of a rigidly rotating perfect fluid with constant mass density (cf. [p]], @). Within 
Einstein's theory of gravitation the problem is even more complicated, mainly 
because of the influence of the body's rotation onto the gravitational field (re- 
lated to the so-called gravitomagnetic potential). Fortunately, for the exterior 
equations a powerful solution technique (the 'inverse scattering method' known 



from soliton theory) exists in the case of axial symmetry and stationarity (|15 - 
[|18|1 ). It does not apply, unfortunately, to the interior equations. However, there 
is an interesting limiting case: For infinitesimally thin disks the interior prob- 
lem 'shrinks' to boundary conditions for the exterior solution, and the solution 
technique mentioned can be utilized for solving the global problem. 

A rigidly rotating disk of dust is the universal limit of rigidly rotating perfect 
fluid configurations as p/e — > (p denotes the pressure and e the mass-energy 
density), cf. the disk limit of the Maclaurin sequence. This disk is interesting for 
two reasons: On the one hand, it represents, in a sense, the simplest model of a 
self-gravitating rotating body (with no interaction except gravitation). On the 
other hand, it may serve as a crude model for astrophysical disks, for example 
galaxies (with the stars considered as dust grains). Of course, normal galaxies 
are sufficiently well described by Newton's theory of gravitation. Nevertheless, 
the relativistic model might be interesting, e.g. in the context of quasars. An 
approximate solution was presented by Bardeen and Wagoner ||, @. 

The rigorous solution of the problem of the rigidly rotating disk of dust ( || - 
0) seems to provide the first example of an exactly solvable rotating-body prob- 
lem within general relativity, apart from the Kerr solution describing a rotating 
black hole. 

The lecture is organized as follows: In Section 2, a brief sketch of the method 
for solving boundary value problems of the axisymmetric, stationary vacuum 
Einstein equations is given. In section 3 the boundary value problem related to 
the dust disk and its solution are discussed. The black hole limit of the solution 
is investigated in some detail. 

I would like to emphasize that most of the material of this lecture is based on 



the joint work with Gemot Neugebauer and Andreas Kleinwachter (H - [12 
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2 Solution of boundary value problems to the 
axisymmetric, stationary vacuum Einstein 
equations 

2.1 The linear system 

The axisymmetric, stationary vacuum Einstein equations (equivalent to the so- 
called Ernst equation, see |13| , [jHfl ) are the integrability condition of a related 
linear system ( ||15|| - ||TH[). Neugebauer's form [H5| of the linear system reads 
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z, A) is a 2 x 2 - matrix function depending on 

2 = P + *C, z = P~K 

and 

A = 



' K — iz 



(1) 
(2) 
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where K is an additional (complex) parameter, called the 'spectral parameter', 
which does not depend on the coordinates p and (. (p and ( are cylindrical 
coordinates, p is the distance to the symmetry - [( -] axis.) A bar denotes 
complex conjugation. The scalar functions M and do not depend on A: 



M = M{z,z), N = N(z,z). 
The integrability condition 



(5) 



(6) 



leads to a first order system^ of nonlinear partial differential equations for M, N, 
M and iV which is equivalent to 



M 



f + f 



N 



f, 



f + f 



(7) 



1 To obtain this system, one has to use the relations 

A 



V-l), A, 2 -=A(A 2 



following from (TJ) . Comparing the coefficients of different powers of A in the resulting equations 
(they must be valid for all Kl) one obtains 

f 



N,g = N(M - N) 
and the complex conjugate relations 



(N + M), M, s = M(N - M) - — (M + N) 
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and the Ernst equation 



W)(U P +f«+- p U) = f, 2 P +f,l (8) 

for the complex function f(p, (), called the Ernst potential. (The relation to the 
metric can be found in section 3.1.) 

The existence of such a linear system with a spectral parameter allows for the 
construction of exact solutions of the corresponding nonlinear partial differential 
equation (here: the Ernst equation), e.g. by means of Backhand transformations 
flT7| ] , |rj| . These solutions contain an arbitrary number of free parameters. More 
importantly, it is even possible to construct solutions containing free functions. 
In this way, one can solve, in principle, initial and/or boundary value problems. 
This method was discovered by Gardner, Greene, Kruskal and Miura in 1967 p(J 
as a method for solving the Cauchy problem of the Korteweg-de Vries (KdV) 
equation. The term 'inverse scattering method' comes from the fact that one 
step of the solution procedure consists in solving an inverse scattering problem 
for the one-dimensional stationary Schrodinger equation which plays the role of 
one part of the linear system related to the KdV equation. 

The general idea behind is the discussion of the matrix function $ as a function 
of the complex spectral parameter K or, in our case, as a function of A. (p and 
( play the role of parameters in this context.) It is possible to obtain solutions 
containing free functions via the solution of related Riemann-Hilbert problems 
in the complex A-plane. This leads to linear integral equations, cf. pi] . 



2.2 The Riemann-Hilbert technique 

A quite general solution $ of the linear system (|J), (fj) can be obtained by solving 
a matrix Riemann-Hilbert problem: This problem consists in finding a $(A) that 
is holomorphic for all values of A in the complex plane except those which lie on 
some closed curves T and T' defined by (£|) and 

k e r K , (9) 

with IV being a closed curve in the complex f^-plane which is symmetric with 
respect to the real axis. (There exist two curves in the A-plane since X(K) is 
double-valued; A 6 T <^ —A G V.) On T and T' the following jump relations 
shall be satisfied: 

<S> l = $ e C(K) onT, $j = $ e C'(K) onT, (10) 

where $j and $ e denote the values of $ that appear by approaching the contour 
from inside and outside, respectively. It is assumed that the jump matrices C(K) 
and C'(K) do not depend on z and z. As a consequence, the expressions $, z 
and $,2 $ _1 do not jump on T and T'. Moreover, together with some additional 
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assumptions in case of zeros of the determinant det$(A), one can show that 
$, z $ _1 and $ _1 are holomorphic functions of A everywhere except at the 
points A = oo and A = 0, respectively. There, in agreement with ([1]) and (Q), 
simple poles occur - provided 

o^r, oo^r. (ii) 

Some constraints on the jump matrices C, C and a suitable normalization con- 
dition ensure that $ _1 and have exactly the structure as in ([[]), @, 
and one can read off the M's and iV's or calculate the Ernst potential f{p,()- 
The solution of a matrix Riemann-Hilbert problem can be found via a system of 
linear integral equations. In this way a solution of the Ernst equation is obtained 
which depends on free functions (some elements of the jump matrices which can 
be choosen arbitralily) . Normally, this solution is regular for all values of p and ( 



satisfying the condition flTT|) , i.e., a curve £ in the p-£-plane has to be excluded 
[cf. (I)]: 

S: p=\SK\, ( = ^K, KeT K . (12) 

This defines the surface of a body of revolution. One can try to solve boundary 
value problems, e.g. of the Dirichlet type, with boundary data given on E. The 
inverse scattering method would then consist of the following steps: 

1. Determination of the jump matrices C(K) and C'(K) from the boundary 
data. [The contour Yk follows from the surface £ according to (p"2|)1. 

2. Solution of the Riemann-Hilbert problem. 

3. Calculation of f(p, Q from $(;z, z, A). 

The first step is the most difficult one. It has to be solved by considering the 
linear system ([!]), (Q) along the boundary S. (In the case of the application of 
the method to the solution of the Cauchy problem of the KdV equation the first 
step is simpler and consists in solving a 'direct' scattering problem: One has to 
determine the 'scattering data' for a given potential.) 

The second step consists in the solution of a system of linear integral equa- 
tions. (In the KdV case it corresponds to the inverse scattering problem: One 
has to reconstruct a potential from scattering data. This leads to the famous 
Gelfand-Levitan-Marchenko integral equation.) 

The third step is trivial and provides us with the desired solution of the 
boundary value problem. 

In the next section we discuss the problem of the rigidly rotating disk of dust, 
which is the first example of a successful application of the procedure outlined. 
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3 The rigidly rotating disk of dust 



3.1 The boundary value problem 

The disk of dust is characterized by the following energy-momentum tensor: 

T ik = eu l u k , u l = e- y (f + Slrf ) , (13) 

with the mass-energy density e and the four-velocity u l . £ l and rf are the Killing 
vectors corresponding to stationarity and axisymmetry, respectively The (pos- 
itive) scalar exp(— V) follows from u l Ui = —1, and Q is the angular velocity as 
measured by an observer at infinity. Rigid rotation means 

Q = constant. (14) 

The line element can be written in the Weyl-Lewis-Papapetrou form 

ds 2 = e~ 2U [e 2 \dp 2 + d( 2 ) + p 2 d V 2 } - e 2U (dt + a d V ) 2 , (15) 

< p < oo, — oo < ( < oo, (16) 

where exp(2U), exp(2/c) and a depend on p and ( only. The Killing vectors, in 
these coordinates, are given by £ l = 51, rf = 5 % . (Note, that we use units where 
the velocity of light c as well as Newton's gravitational constant G are equal to 
1.) The disk is defined by ( = 0, p < po; po being the (coordinate) radius of the 



infinity: 
/-I 



disk (£): 



Pl — e 2V 



P 



Po 



Figure 1: The boundary value problem, 
disk. The mass density e(p, Q can be written formally as 

6 = 5(()e u - k a p (p), 



(17) 
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where S(Q is Dirac's delta function and u p (p) is the (proper) surface mass-density. 
The boundary conditions can easily be obtained by integrating the Einstein equa- 
tions across the disk, taking into account that the metric coefficients are continu- 
ous, but their normal derivatives may jump. Together with the symmetry of the 
problem with respect to the plane ( = one obtains 

e 2U' = e 2V = constant ^ a ' )f = 0j k', c = as C = 0, p < po. (18) 

The primed quantities refer to the corotating frame of reference defined by 

p' = p, (' = (, cp' = <p - fit, t' = t. (19) 

The constant Vq is related to the relative redshift zq measured by an observer at 
infinity for photons coming from the center of the disk: 

z = e~ Vo - 1. (20) 

[Note that the V in (|13|) turns out to be constant and equal to Vq.] The surface 
mass-density can be calculated (after one has solved the global problem) by 



(21) 

C=o+ 



To formulate the global problem as a boundary value problem for the Ernst 
equation, we need the connection of the metric functions with the complex Ernst 
potential f(p, (): 

f = e 2U + ib, (22) 

i.e., the real part of the Ernst potential is just equal to exp(2U). The imaginary 
part b is related to the gravitomagnetic potential a: 

a, p = pe~ 4U b, ( , a, c = -pe~ 4U b, p . (23) 

Hence a(p, () can be calculated via a line integral from the Ernst potential and 
its derivatives. This integral is path-independent. [The integrability condition 
is satisfied as a consequence of (§).] Similarly, the metric function k(p, () follows 
from 

k, p = p[U, 2 p -U* +\e-* u (b, 2 p -b, 2 c )], k, c = 2p(U, p U, c +\e-* u b, p b, c ). (24) 

Note that, for regularity reasons, a = and k = on the symmetry axis (p = 0), 
and we obtain, for example, 

p 

pe- 4U b <( ip, (25) 

o 
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k = j p[U,j -U, a c +\e-* u (b* -b* )]d~p. (26) 
o 

[In the integrands, one has U — U(p, £) and b = b(p, ().] 

Now, from (|1"8"D, (^) and fl23|) we conclude that the boundary condition for 
the 'corotating' Ernst potential /' is simply given by [fj 

/' = e 2V ° = constant as C = 0, p < po- (27) 

It should be noted that we have combined the result b', p — and the freedom of 
adding an imaginary constant to the Ernst potential, to set b' = in the disk. 
The simple condition ( p?[) for /' in the corotating system corresponds to a quite 
complicated, nonlocal boundary condition! for / in the original system. On the 
other hand, the condition of asymptotic flatness is much simpler for /: 

/ -»• 1 as p 2 + C 2 -> oo. (28) 

This condition ensures that U — > 0, k — > and a — > at infinity, i.e., the line 



element (|T5 ) becomes Minkowskian (in cylindrical coordinates). 

The global problem one has to solve, is to find a (or: the) solution of the 
Ernst equation (Q) which satisfies (|2"TD and fl28|), and which is regular everywhere 
outside the disk, cf. Figure 1. 

The solution depends on two parameters only. One can choose, e.g., Vq and 
Po, or Vq and Q. A relation Q = Q(Vo,po), see next section, follows from the 
regularity condition at the rim of the disk. 

3.2 The solution 

The solution was obtained by applying the method outlined in Section 2. The 
curve Tk in the complex i^-plane is a part of the imaginary axis (from — ip to 
+ipo) corresponding to the surface E, i.e. the disk, cf. (0) and Figure 1. 

The first step of the solution procedure (the determination of the jump ma- 
trices) lead to the 'small' integral equation || which could be solved in terms of 
elliptic functions |J. The 'big' integral equation corresponding to the second 
step (the solution of the matrix Riemann-Hilbert problem) could be solved in 

2 The coordinate transformation ( |l9| ) preserves the form of the line element (|l^). One obtains 



z 2U ' = e 2U [(l + naf - n 2 p 2 e~ iU ], (1 - n a ')e 2U ' = (1 + n a )e 2U ', e 2k '- 2U ' = e 2k - 2U . 



To get b' , one has to integrate the (primed) relations (p3[). That means, the boundary condition 
for / is given by complicated, nonlinear and nonlocal relations corresponding to /' = e 2Vo . Note 
that these boundary conditions are equivalent to those applied by Bardeen and Wagoner S, 
PJ. Vq is identical with the parameter v c of Bardeen and Wagoner. 
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terms of hyperelliptic functions and lead to the following result for the Ernst 
potential 0: 



/ = exp 



l K 2 dK 



K, 



K 2 dK 



V 2 



Ki 



with 



Z = y/(K + iz) (K — iz) (K 



K\){K 2 



(3Mfi < 0), K 2 



-KI), 



The real (positive) parameter p is given by 

V = 2n 2 p 2 e- 2Vo . 

The upper integration limits K a and Kb in ( p9|) have to be calculated from 



k„ 



A'l 



dK 
IT 



K h 



Ki 



dK 



KdK 



K h 



+ 



KdK 



Vl, 



(29) 

(30) 
(31) 

(32) 
(33) 



Kx 



Ki 



where the functions Vq, V\ and v 2 in (|33"D and (|29| ) are given by 



+ipo 



v 



Z~i 



+ipo 



dK, vi 



H 



+ipO 



KdK, v 2 



H 



K 2 dK, 



H 



/iln J I + p?{\ + K 2 /plY + p(l + K 2 /pl 



m.H = o), 

mp 2 ^l + p 2 (l + K 2 /p 2 ) 2 

(JftZi < for p, ( outside the disk). 



(34) 

(35) 
(36) 



Z x = y(K + iz)(K — iz) 

In (|33p one has to integrate along the imaginary axis. The integrations from 
Ki to K a and K 2 to K b in (|29|) and (|33D have to be performed along the same 
paths in the two-sheeted Riemann surface associated with Z(K). The problem 
of finding K a and K^ from (|3^) is a special case of Jacobi's inversion problem. 
It generalizes the definition of elliptic functions and can be solved in terms of 
hyperelliptic theta functions ([fm, [0' see a ^ so (E3 HH)- Using a formula 
for Abelian integrals of the third kind derived by Riemann (see [[24]]) it is also 
possible to express the Ernst potential / directly in terms of theta functions [|T2]|. 
On the symmetry axis (p = 0) and in the plane of the disk (£ = 0) all integrals 
in (^) and (|33[) are reduced to elliptic ones f|. 

The solution (|29| ) satisfies the boundary conditions ( p7|) and (p8|), has a posi- 
tive surface mass- (particle number-) density (vanishing at the rim of the disk), 
and it is regular everywhere outside the disk - provided 



< ji < p = 4.62966184. . . 



(37) 
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(for p > p one or more singular rings appear in the plane ( = 0, outside the 
disk). The interesting behaviour for p — > p will be discussed in the next section. 

Note that the solution in the form (^) - ( 136|) depends on the parameters po 
and p only. Since exp(2U') = exp(2U) on the symmetry axis (p = 0), one can 
calculate the parameter Vo [cf. (|T8|)1 from 9?/(p = 0, £ = + ). The result is ||: 

Vb = -^sinh' 1 |p + 4 — rx-r; , — ^77^ 2-}> (38) 

2 [ p[/(p);fp 2 -4,fp(l + p 2 /9)]-fpJ 



/(,) = i/i^i^£ (39) 



(p is the Weierstrafi function ^), i.e., Vo depends on p alone. The range < 
p < Po corresponds to > Vo > — oo. In this range, the relation (|38|) can be 
inverted uniquely to give p(Vo). [po is the first zero of the denominator in (p8|) .1 
Then, from the definition ( |32] ) one obtains the relation f2(Vo,po). (Without loss 
of generality, we assumed Q > 0; the solution for negative Q is simply given 
by /.) Alternatively, one can use Vo and Q as the primary parameters, with 
Po = Po(Vo,^). 

3.3 Discussion: 

^From the Newtonian limit to the black hole limit 

For < p < po, the solution ( p9f) can be expanded in terms of p 1 / 2 : 

oo 

/ = 1 + J2 fn /i ( " +1)/2 = 1 + fi A* + CV /2 )- (40) 

n=l 

The /„ (n = 1,2, . . . , oo) are elementary functions of p and £. (This series 
corresponds to the Bardeen- Wagoner expansion^ .) The Newtonian limit (//< 1) 
is represented by f\: 

fx = -\ { \ co^ 1 £ + K - (£ 2 + |) cot- 1 £] (1 - 3p 2 ) } , (41) 
with elliptic coordinates £ and 77: 



p = Po V / l + ev /1 -^ C = Po£? (0<£<oo, — 1 < 77 < 1) (42) 



3 Thc Weierstrafi function #2?#3) is defined by 

00 

dt 



s v 7 ^ 3 - git - 53 

4 The expansion parameter 7 used by Bardeen and Wagoner ||, Q is related to by 7 
]_ _ e VoW = M / 2 + 0( M 2 ). 
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(£ = is the disk). The Newtonian potential Un is obtained from 

or j i 

9u = -»/ = -(1 + A ai) + 0(^ 2 ) = -(1 + + O(-), (43) 

i.e., 

U N = n 2 p 2 f 1 . (44) 

This is exactly the Maclaurin solution. [Note that we have reintroduced the 
velocity of light c into Eq. fl4~3~D . From ( j32|) and (|3~8"D one obtains /i = 2VL 2 p 2 ] /c 2 + 
CV).] 

With increasing /i, characteristic deviations from the Newtonian solution oc- 
cur. This concerns, e.g., the radial distribution of the surface mass-density 0,0. 
The gravitomagnetic potential leads to dragging effects and, for // > 1.68849 . . ., 
even to the formation of an ergoregion |3J], [|TTJ. Some illustrations can also be 



found in [12 



However, the most striking difference to the Newtonian case is the following 
i, & 

For given angular momentum the mass of the disk is bounded by 

M 2 , N 

— < 1, (45) 

where M and J denote the total (gravitational) mass and the ((^-component of 
the) total angular momentum, respectively (c = G = 1 again). The equality 
M 2 j J = 1 is reached in the limit \i — > the black hole limit. 

For ji — > /i , one has V — > — oo, cf. fl3"E|). As a consequence of (§^), for 
nonvanishing Q, this results in 

Po -> 0. (46) 
For p 2 + ( 2 ^ 0, we obtain from ( p9|) in the limit /i — > yU 

2fir — 1 — i cos i? , . 

' = 2flr + i-«™« < r > »>• < 47 > 

p = rsin , #, £ = rcosi9 (0 < $ < 7r). (48) 
This is exactly the (r > part of the) extreme Kerr solution ^ with 

M = 2S' J -w-' (49) 



To derive (|47|) from (J29|) , let us first rewrite ([29|) and ( p3| ) in the equivalent form 

K a \ K a K a 

, f K 2 dK f dK f KdK 

7=exp<( / — v 2 >, ~ = V °' / — Z~ = Vu 

VK b ) K b K b 

with 

f K n dK , 
v n =v n - / — - — (n = 0,1,2). 
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i.e., 

M 2 , . 

~J = ^ 
(Remember that we assumed Q > 0.) In the coordinates used, the horizon of the 
extreme Kerr black hole is just given by (the excluded) r = 0. Q plays the role of 
the 'angular velocity of the horizon'. In the corotating system (|19|), one obtains 



2(1 + z cos??) 2 „ ' 

+ sin 2 ■& 



2Vtr + 1 — i cos $ 



(51) 



It can be seen, that the boundary condition (|27f) with V — > — oo is indeed satisfied 
on the horizon (r = 0). 

A completely different limit of the space-time, for \i — > /i , is obtained for 
finite values of r/p (corresponding just to the previously excluded r = 0). There- 
fore, we consider a coordinate transformation Q 

r = re- Vo , <p = (p-ttt, = t = te Vo . (52) 

(Note that finite r/p correspond to finite f in the limit.) For fi < /i , this is 
nothing but the transformation to the corotating system ( |19|) combined with a 
rescaling of r and t. The transformed Ernst potential / is related to /' according 
to/ = /'exp(-2V ),i.e., 

/ /' 

— = — as n < fj, . (53) 

However, for \i — > //q, the solutions /' (finite r) and / (finite f) separate from 
each other. (A similar phenomenon has been observed by Breitenlohner et al. for 
some limit solutions of the static, spherically symmetric Einstein- Yang-Mills- 
Higgs equations f2"?fl .) For finite r, the extreme Kerr solution ([51]) arises, while 
finite f lead to a solution which still describes a disk. This solution (which can be 
expressed in terms of theta functions) is regular everywhere outside the disk, but 
it is not asymptotically flat, i.e., it can be considered as a cosmological solution. 
The space-time structure of both solutions (/' and /) coincides at r = (the 
horizon) and r —>■ oo (spatial infinity). The relation (Q) survives in the form 

/ /' 

lim — = lim — as fi — > /io- (54) 

twoo f i — >0 T 



{Kb is now on the other sheet of the Riemann surface.) In the limit fj, — > fio one obtains for 



r > 0, using (38) 



2£! 7ricos?9 iri 

V0 = 7T1— » v l=-7Ti V 2 ={) 

r 2r z Zr 

(modulo periods). In the above integrals from Kb to K a , because of K\ — » 0, if 2 — ► [cf. (|31_ 
Z can be replaced by Z = K 2 \J (K + iz)(K — it). Hence, all integrals become elementary and 
the result p7|) can easily (and uniquely) be obtained. 
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Accordingly, for fi — > //o and r — > oo, 

"2(1 + icos^) 2 



2-2 



/ - /„ = -0 2 f 



1 — Z COS X? 



+ sin 2 1? 



(55) 



Note that / os belongs to the family of solutions to the Ernst equation of the 
type / = r k Yf t (cos , d) presented by Ernst PSJ. The corresponding asymptotic line 



element is given by the following exact solution of the vacuum Einstein equations: 

ds 2 = e^[e 2i (rfr~ 2 + r~W) + f 2 sin 2 t? # 2 ] - e 2& (dt + a d<p) 2 , (56) 

ofr ^9 9 cos 4 $ + 6 cos 2 ^9 — 3 . , 

e 2U = Q 2~2 ^ ^ ( 5? ) 

COS 2 1? + 1 

- 2 C0S2 *" 1 (58) 



fi 2 r cos 4 t? + 6 cos 2 - 3 ' 

e 2 ^ = I( cos 4 ^ + 6cos 2 ^-3). (59) 

These analytical results prove the conjectures formulated by Bardeen and Wag- 
oner H| on the basis of their numerical results. 

Let me conclude with a quotation from Bardeen and Wagoner page 411): 
T/ie picture we have developed of the extreme relativistic limit of a rotating disk 
is that it becomes buried in the horizon of the J = M 2 Kerr metric, surrounded 
by its own infinite, non asymptotically flat universe. As approached from the 
r > Kerr region, the disk represents a "singularity" in the horizon, since the 
whole range < f < oo, over which there exist considerable changes in the local 
geometry, corresponds to an infinitesimal range of affine parameter for a typical 
photon which reaches the horizon from the outside. ' 
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